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We study number-phase uncertainty in a laser-driven, effectively four-level atomic system under 
electromagnetically induced transparency (EIT) and coherent population trapping (CPT). Uncer- 
tainty is described using (entropic) knowledge of the two complementary variables, namely, number 
and phase, where knowledge is defined as the relative entropy with respect to a uniform distribu- 
tion. In the regime where the coupling and probe lasers are approximately of equal strength, and 
the atom exists in a CPT state, there is coherence between the ground states, and correspondingly 
large phase knowledge and lower number knowledge. The situation is the opposite in the case where 
coupling laser is much stronger and the atom exists in an EIT state. We study these effects also in 
the presence of a higher-order nonlinear absorption, which is seen to produce a dephasing effect. 



which holds accurately when the number n is large. When the number uncertainty is small, A<j) should be large 
but the limit on <f> is 2tt and the uncertainty relation breaks down. Also, it is not possible to introduce a hermitian 
phase operator because n is bounded below [T]. One way to get around the difficulty is to introduce a probability 
distribution V(4>) of phase, where the operator corresponding to phase is not a conventional projective measurement, 
but a positive operator- valued measure (POVM) in particular, given by projectors to the atomic coherent states 
\0,(p). As the expectation value of the dipole moment operator for these states is j(sm8)e^ 1 ^ , where j is the total 
angular momentum of the electron, measurement of this POVM yields phase information. Phase 4> is an important 
quantity in atomic coherences and in the interferometry based on such coherences [3l |4] . 

In this paper, number-phase complementarity in an atomic system interacting with a laser field is studied via an 
information theoretic treatment. A particular context in which this could be studied is in spin squeezed atomic 
ensembles, which can be created either by transfer of squeezing from squeezed light fields or by the presence of 
nonlinearities |S]. Here we consider an atom (e.g., Rb) with a four-level electronic hyperfine manifold under the action 
of two laser fields. This system exhibits several interesting nonlinear phenomena such as electromagnetically induced 
transparency (EIT), coherent population transfer (CPT) and three-photon absorption (TPA), prompting us to study 
the number-phase complementarity in these regimes. 

The plan of the paper is: Section [TT| describes the atom laser sytem. Section [TTl] briefly introduces the concept of 
phase distribution; in Section |IV| the information theoretic formulation of complementarity is presented in terms of 
an upper bound to a (convex) sum of knowledge of the complementary observables, number and phase. In Section 
|Vj the entropic uncertainty relations are used to study number-phase complementarity in EIT and CPT systems. In 
Section |VB| the effect of higher order nonlinearity on the complementaristic behavior is studied. In Section |VI[ we 
make our conclusions. 
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The uncertainty relation in number and phase is 



AnAcj) > 1/2, 
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FIG. 1: The atom-laser system: couplings provided by the lasers. The coupling laser gives rise to two Rabi frequency couplings 
(3 — 3' and 3 — 4') and the probe to one (2 — 3'). The coupling that does not participate in EIT formation (3 — 4') gives rise to 
third order nonlinearity in probe absorption. There is nonlinear absorption in the coupling f2 34 / also. This coupling is called 
the signal beam in the XPM scheme. 

II. BRIEF DISCUSSION OF THE ATOM LASER SYSTEM 

In this study we investigate the complementarity in atomic states driven by unsqueezed light. The isotope under 
study, 85 Rb, has a nuclear spin of 5/2. The electronic levels of interest in conventional atomic spectroscopy notation, 
are the ground hyperfine levels 5Si/ 2 F — 2,3 and the excited states 5P 3 / 2 F = 3', 4'. The selection rules permit 
electric dipole transitions that satisfy AL = ±1 and AF=0, ±1. Two lasers, coupling and probe, at frequency 
loc and up, detuned by Scoupimg (Sc) and 5p ro be (dp) from the two transitions 5Si/ 2 F = 3 — > 5P 3 /2-F = 4' and 
5Si/ 2 F = 2 — > 5P 3 / 2 P — 3' with Rabi frequencies H 33 /, Sl 34 / and ^23', are incident on the atom, see Fig. (ITJ). 

A three-level A system is formed by the levels 5Si/ 2 F — 2,3 and 5P 3 / 2 F = 3'. The electron may make a transition 
between the two ground hyperfine states via the excited level F — 3'. However, depending on the Rabi frequencies 
of the two lasers, either a cancellation of transition amplitudes leading to EIT or the formation of a superposition 
dark state leading to CPT, is expected when the detunings of the lasers satisfy Sc = 5p. Due to the Rabi frequency 
coupling SI34' , there is a three photon absorption taking place and the probe laser sees a nonlinearity [3 IE]- There 
is nonlinear absorption in the signal beam also @ [3] . 

In [10] , the four level density matrix equations in the presence of the two beams and three Rabi frequency couplings 
were solved. There the Lindblad terms corresponding to spontaneous decay were used. 

To make a study of the four level system in number-phase variables we recast the relevant 4 levels of the atom in 
terms of a pseudospin of spin 3/2. The projections of the spin can take any of the four values s = —3/2, —1/2, 1/2, 3/2, 
which we take to represent the four electronic levels, and are mapped to F = 2, 3, 3', 4', respectively. It can be shown 
that a different assignment merely constitutes a re-labelling of vectors, and does not qualitatively alter the conclusions 
that follow. 

Following the notation of (llj we evaluate the phase distribution function P(0), from which we obtain phase 
knowledge R[0] and number knowledge R[m], defined below. The level populations are reflected in the angular 
momentum basis (number, m) while the coherence between levels is reflected in the phase variable (j). 
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III. PHASE DISTRIBUTION 

The quantum description of phases [T2J [13] has a long history [1] IT4HT7] . Pegg and Barnett [16] , following Dirac 
[14j . carried out a polar decomposition of the annihilation operator and defined a hermitian phase operator in a 
finite-dimensional Hilbert space. In their scheme, the expectation value of a function of the phase operator is first 
carried out in a finite-dimensional Hilbert space, and then the dimension is taken to the limit of infinity. However, it 
is not possible to interpret this expectation value as that of a function of a hermitian phase operator in an infinite- 
dimensional Hilbert space [2[ 118) . To circumvent this problem, the concept of phase distribution for the quantum 
phase has been introduced [T51 \W\ . In this scheme, one associates a phase distribution to a given state such that the 
average of a function of the phase operator in the state, computed with the phase distribution, reproduces the results 
of Pegg and Barnett. 

For an atomic system, the phase distribution 7-* (</>), 4> being related to the phase of the dipole moment of the system, 
is given by [3J 

V{<f) = 2] £± £ d9 sm(8)Q(0, 0), (2) 

where V((j)) > and is normalized to unity, i.e., d<j)V((j)) = 1. In the above, j is the pseudo angular momentum 
of the atom. Here Q(0, 4>) is defined as 

q(M) = <Mp"M>, (3) 

where \0, cf>) are the atomic coherent states [501 [2T] defined by Eq. Q in terms of Wigner-Dicke states [25], which are 
the simultaneous eigenstates of the angular momentum operators iF and J z ; 



E 



2j 
j + m 



(sm(6/2)y +m (cos(6/2)y- m \j, m)e^ +m ^. (4) 



It can be shown that the angular momentum operators Jj, and (obtained by rotating the operators J x , J v and 
J z through an angle 6 about an axis n = (sin <f>, — cos </>, 0)), being mutually non-commuting, obey an uncertainty 
relationship of the type > \{J^). Atomic coherent states (Eq. 4| are precisely those states that saturate 

this bound, hence the name, in analogy with radiation fields [3]. For two level systems, they exhaust all pure states, 
whereas for larger dimensions, this is no longer true. Using Eq. ^ in Eq. ([2]), with insertions of partitions of unity 
in terms of the Wigner-Dicke states, we can write the phase distribution function as [11] 



4-7T 



V(^) = T dOsine Y, y.-ii.n)U.n\int)\j.n,)(j.m\».n). ( ) 



n,m——j 



The phase distribution "P (</>), taking into account the environmental effects, has been studied in detail for QND as well 
as dissipative systems in [TT] [23] for physically interesting initial conditions of the system S, i.e., (a) Wigner-Dicke 
state, (b) atomic coherent state and (c) atomic squeezed state. 

In our mapping scheme, Wigner-Dicke or excitation states are thought of as 'number states', thereby making J z 
the 'number observable', whose distribution p(rn), given below as 

p{m) = (j, m \p s (t)\j,m), (6) 

is considered as complementary to V((f>) [3J- 

IV. INFORMATION THEORETIC REPRESENTATION OF COMPLEMENTARITY 

Two observables A and B of a <i-level system are called complementary in quantum mechanics if measurement of 
A disturbs B, and vice versa [211 ESj- Complementarity is related to the Heisenberg uncertainty principle, which says 
that for any state "0, the probability distributions obtained by measuring A and B cannot both be simultaneously 
peaked if A and B are non-commuting. Heisenberg uncertainty is traditionally expressed by the relation 

A^AA^B> ±\([A,B])i,\, (7) 
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where (IS^A) 2 = (A 2 )^ — ((A)^) 2 . However, this representation of the Heisenberg uncertainty relation has the 
disadvantage that the right hand side of Eq. ^ is not a fixed lower bound but is state dependent. Further, the form 
of Eq. ([7| is not invariant when A or B is scaled by some numerical factor, though physically we would want the 
measure of non-commutativity to be scale- invariant. 

The information theoretic (or "entropic") version of the Heisenberg uncertainty relationship [24 -26 , which uses 
Shannon entropy of measurement outcomes, instead of variance, as a measure of uncertainty |27[ 128]. overcomes both 
these problems. 

The relative entropy associated with a discrete distribution f{j) with respect to a distribution g(j) defined over the 
same index set, is given by 

*C« = £/C*)i°g(^)- (8) 

S(f\\g) > can be thought of as a measure of 'distance' of distribution / from distribution g , where the equality 
holds if and only if f(j) — g(j) [27 . Consider a random variable F with probability distribution /. We will define 
R(F) as the relative entropy of / with respect to the uniform distribution which is 4, for a system of dimension d (as 
the system has equal probability of being in any of the d states), i.e., 

R(F) = R[f(j)} = ]T f(j) log(d/(j)). (9) 

3 

As a measure of distance from a uniform distribution, which has maximal entropy, R(F) can be interpreted as a 
measure of knowledge, as against uncertainty, of the random variable described by distribution /. Following [29], we 
recast the Heisenberg uncertainty principle in terms of relative entropy as 

R{A) + R(B) < logd, (10) 

where d is the (finite) dimension of the system. The Hermitian observables A and B are said to correspond to mutually 
unbiased bases (MUB-s) if any eigenstate of one of the observable can be written as an equal amplitude superposition 



of all the eigenstates of the other observable [25]. Physically, Eq. (10) expresses the fact that simultaneous knowledge 
of A and B is bounded above by logd, and that the probability distributions obtained by measuring A and B on 
several identical copies of a given state cannot both peak simultaneously. 
Eq. QSJ) has a natural extension to the continuous case, given by 

S(f\\g) = J d pf(p)log(^y (11) 

As in the discrete case, we define R(f) as relative entropy setting g(p) to a continuous constant function. In particular, 
the relative entropy of V((f>) with respect to a uniform distribution ^ [TTJ [23], corresponding to the case where phase 
is completely randomized, over <p is given by the functional 

R[P (</>)] = f #P(^)log[27rP(^)], (12) 
Jo 

where the log(-) refers to the binary base. 

For the phase POVM measure <j>, because of the non-orthogonality of the states \9,<j>), the concept of an eigenstate of 
an observable is weakened to that of a minimum uncertainty state, which corresponds to maximum phase knowledge. 
On the other hand, number is a regular, Hermitian observable. This leads to the concept of MUB being replaced by 
that of a quasi-MUB. 

Two variables A and B form a quasi-MUB if the minimum uncertainty state of A is a maximum uncertainty state 
of B, but the minimum uncertainty states of B are not necessarily maximally uncertain in A. The number-phase 
complementarity in our four- level system indeed shows such quasi-MUB character [25]. This can be seen by noting 
that for the Wigner-Dicke states \j, to) , the phase distribution is pj] 

= 2J ^( j % 2 ) B \j + to + l,j - rn + 1] = ±, (13) 



where B stands for the Beta function. Thus, it follows via Eq. (12 1 that the knowledge R^ vanishes. On the other 
hand, it can be shown that in the present case the states which minimize R^ are the not Wigner-Dicke states. However, 
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it can be shown that they are number states in the d = 2 case. To see this, we observe that if V(4>) is constant, then 
in Eq. ([H]), each term in the summation, which is proportional to e l ( m -™)0 j must individually be independent of <j>. 
Since 4> is arbitrary, this is possible only if m = n, i.e., the state p s is diagonal in the Wigner-Dicke basis. 

Thus J z and <fi form a quasi-MUB because of the POVM nature of </> [27]. For a POVM, knowledge R of the minimum 
uncertainty state need not be log(<f) bits, essentially because outcomes are non-sharp owing to non-orthogonality |30] 
of the corresponding measurement operators [29]. Hence, the plain summation over knowledges in the entropic 
formulation of the uncertainty relation is replaced by a suitable convex sum. 

Thus, for a POVM in the four level case such as phase in our case, we have that for minimum uncertainty states 
R < 2 bits. To compensate for smaller values of R^ and to be able to define coherent states in a four level system 
that are not Wigner-Dicke states, we introduce the parameter [i2 > 1, to obtain inequality: 

R s (n 2 )=»2R<i> + Rm<2 (14) 

over all states in C 4 . Our strategy is to numerically search over all states in this space, other than the Wiger-Dickc 



states, where R m = 2 and R^ — 0, saturating the inequality (14 1 trivially- in order to determine the largest value of 



fj,2 such that this inequality is just satisfied. By this method, we find /12 = 1.973. 



An alternative formulation of Eq. (14), which we do not use here, is to rewrite the summation in the right hand 
side as a convex sum: 

+ (1 - p)R m < c, (15) 

where p = /i2(l + M2) an d c = 2(1 + P2) 1 , which can now be considered as a generalized uncertainty relation in 
the sense of [31]. For the present case, p = 0.66 and c = 0.67. 



V. NUMBER-PHASE COMPLEMENTARITY IN CPT AND EIT SYSTEMS 



By coherences, we mean the off-diagonal elements of p in the number representation, while by populations we mean 
the diagonal elements. Because, as we saw, a number state has a flat distribution of phase, so does a mixture of 
number states, which is represented by a purely diagonal density operator. Thus, non-vanishing phase knowledge 
implies that there is coherence in the state. 

In this and the following section, we apply the ideas and tools of the preceding sections to the atomic system 
described in Section [IT] The density matrix calculations and the evaluation of the entropies was carried out for a wide 
range of values of parameters 5c, Sp, ^33', ^23', ^34'- 



A. The coupling-probe A-system 

Initially we discuss the case of CPT (both lasers of equal strengths) and then EIT (lasers of very dissimilar strengths). 
To begin with, we omit the coupling $734/ and set coupling f2 2 3' = ^33' = 5. The phase distribution V{(t>) for the case 
of probe and coupling being on resonance with their respective transitions, is non-uniform and peaks at (f> = tt. Thus 
phase knowledge is non-vanishing, as expected. On the other hand, with the coupling laser at resonance, but with 
the probe detuned away from resonance (no CPT or coherence in /O23), the number knowledge remains roughly the 
same but the phase distribution tends to become uniform, suggestive of dephasing. By contrast, making the coupling 
laser much stronger at resonance produces a uniform phase distribution by creation of a number state under EIT 
conditions. These ideas are systematically presented in Figure [2] 

Figure [2] is an array plot depicting the degree of mixedness (1 — Trp 2 ) (left column), phase knowledge (center 
column) and number knowledge (right column) as we proceed from CPT (top row) to EIT (bottom row). Each plot 
has detuning of the probe laser on the abscissa. 

In each regime (row), we find that the system is in a pure state on resonance, and thus amenable to an interpretation 
in terms of complementarity. The following discussion pertains to the system under resonant condition. Comparing 
the different regimes, we find that coherence (phase knowledge) is largest in CPT and number knowledge the least, 
while the opposite is true in the EIT case. Further, the number knowledge being 2 bits in the latter case implies 
that the atom exists in a definite number state. This state can be readily identified with level F — 2, which is 'dark' 
under the stronger 3 — 3' laser. On the other hand, the number knowledge being 1 bit in the former case, with phase 
knowledge being large, implies that the atom exists in an equal weight state ^4= (| 2) + e l *|3)). 

Going off resonance, the situation changes considerably in the CPT case in that phase knowledge vanishes and 
mixedness becomes maximal, implying a transition from a coherent superposition to a statistical mixture. Note 
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FIG. 2: The parameters are f2 33 / (coupling) = S7 23 / (probe) = 5 MHz for the first row, corresponding to a CPT state. For the 
second and third rows, the probe is 2.5 and 0.5 MHz, while coupling beam strength remains the same, corresponding to an 
intermediate and an EIT state, respectively. In all these figures, the a>axis represents ramping of the probe by means of detuning 
from level 3'. The CPT position is at detuning Sp— 1MHz of the Probe laser. On resonance, the figures show complementary 
behavior between number and phase. Going off resonance, phase knowledge is affected, but not number knowledge, implying a 
dephasing effect. This also explains why there is loss of purity in the CPT and intermediate case (first and second rows, where 
there is non- vanishing phase knowledge) but not in the EIT state (third row), where there is complete number knowledge (2 
bits) and no phase knowledge. 



however that the number knowledge remains the same, which signifies that the noise is of a dephasing kind, i.e., the 
phase <j) gets randomized while populations remain unchanged. 

Thus, in going down the central column of the figure, the mechanism of loss of phase knowledge on-resonance is 
the CPT-to-EIT transition, whereas off-resonance it is dephasing when the coupling and probe laser strengths are 
comparable. In the purely EIT case (bottom row), the off-resonance purity is attributed to optical pumping into the 
level 2. 



B. Effect of higher order nonlinear ity on phase knowledge 

Thus far we had considered an atom under the action of two light fields, ^2,3' and ^33'. However, the atom under 
consideration is a multi- level system with additional closely spaced excited levels- 4' and 2'. The transition 3 — 4' is 
an order of magnitude more probable than 3 — 2'. We therefore introduce a term in Sl^i representing the off-resonant 
coupling 3 — 4' of the coupling laser to a nearby excited line. This results in a higher order nonlinearity, the effects 
of which are discussed below. 

In the CPT case, the coupling disrupts the dark coherent superposition state formed during CPT by bringing 
about induced transitions to the state 4'. This lowers the phase information as strength of the signal beam is increased 
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FIG. 3: The effect of the off-resonant coupling 3 — 4' (signal beam) on the CPT state of Figure [2] The coupling and probe 
beams, for all three rows, are taken at strengths 5 MHz. The signal beam strength is taken to be 0.5, 2.5 and 5 MHz respectively 
for the top, middle and bottom rows, respectively. 



( second column of Fig. [3]) as where there is a possibility of spontaneous emission. This is equivalent to measuring 
the CPT state in the basis {\F = 2), \ F = 3)} with a probability determined by ^34/, and thus equivalent to an 
application of a phase damping channel of strength determined by ^34' . This has the effect of randomizing the phase 
in the CPT state, thereby decreasing phase knowledge but not affecting the number knowledge (third column of Fig. 



VI. RESULTS AND CONCLUSIONS 



EIT and CPT states present two contrasting nonlinear phenomena in optics that illustrate the complementary 
behavior of number and phase in atomic systems. A conventional description of this complementarity, based on the 
non-commutativity of these two variables is not possible, as lower-boundcdness prevents the possibility of phase as a 
Hermitian observable. One way out is to represent number by a continuous- variable POVM, described by a probability 
distribution for a given state. Complementarity can then be quantified, among other ways, by expressing the spread 
in the respective distributions by the entropy generated by measurement. The issue of employing the discrete- valued 
Pegg-Barnett phase operator instead of the POVM used here will be discussed elsewhere. 

We have used entropic knowledge, rather than variances, to describe uncertainty in n and <f>. We find that in CPT, 
the coherence between the ground states participating in the dark state is reflected in large phase knowledge and 
about 1 bit of number knowledge. In EIT, where the dark state is a number state, number knowledge is maximal (2 
bits), while phase knowledge vanishes. Thus on-resonance, we see a clear manifestation of number-phase uncertainty. 
Off-resonance, there is a general reduction in coherence, and, in the CPT case, a reduction in phase knowledge due 
to dephasing. A similar phase damping effect is seen also when a higher order nonlinearity is introduced by allowing 



for 3 — 4' transitions. 
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